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Quark gap equation within the analytic approach to QCD 
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The compatibility between the QCD analytic invariant charge and chiral symmetry breaking is examined 
in detail. The coupling in question incorporates asymptotic freedom and infrared enhancement into a single 
expression, and contains only one adjustable parameter with dimension of mass. When inserted into the standard 
form of the quark gap-equation it gives rise to solutions displaying singular confining behavior at the origin. 
By relating these solutions to the pion decay constant, a rough estimate of about 880 Me V is obtained for the 
aforementioned mass-scale. 
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As has been advocated in recent years by a 
number of authors PP, a possible way for gain- 
ing further insight on the infrared sector of QCD 
is to resort to arguments based on the analyticity 
of the underlying dynamics, as captured by the 
corresponding dispersion relations. Jointly with 
the renormalization group formalism, these rela- 
tions constitute the essential ingredients of the so- 
called "analytic approach" [Q to Quantum Field 
Theory (some applications of this method can be 
found in Refs. |ll2M4| b 

One of the most celebrated applications of this 
approach to QCD is related to the running cou- 
pling. In particular, one is able to extrapolate 
the behavior of the strong coupling from the ul- 
traviolet region towards the infrared regime, by 
imposing requirements of analyticity on it or on 
its derivative |HHj . Specifically, if we apply such 
ideas on the perturbative expansion for the one- 
loop (3 function, the solution for the analytic in- 
variant charge that emerges assumes the form [5] 
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The invariant charge has the correct analytic 
properties in the k 2 variable, i.e., unlike its per- 
turbative counterpart, it displays no Landau pole. 
In addition, it has no adjustable parameters other 
than a mass-scale, to be denoted by A, emerg- 
ing, as usual, through the standard mechanism of 
dimensional transmutation. Moreover, the cou- 



pling of incorporates asymptotic freedom and 
infrared enhancement into a single expression (see 
Figure^). Regarding this last point, the invariant 
charge Q has been shown to generate the confin- 
ing static quark- antiquark potential with a quasi- 
linear raising behavior at large distances; specifi- 
cally V(r) ~ rj lnr when r — > oo (r is the dimen- 
sionless distance between quark and antiquark). 
Further appealing features of the analytic running 
coupling Q and its applications can be found in 
Ref. 3 . Evidently, it would be interesting to fur- 
ther scrutinize the advantages and possible limi- 
tations of the analytic charge in the region where 
it aspires to make a difference, namely the in- 
frared sector of QCD. 

Chiral symmetry breaking (CSB) and dy- 
namical mass generation are inherently non- 
perturbative effects, whose study in the contin- 
uum leads almost invariably to a treatment based 
on the Schwinger-Dyson (SD) equations of the 
theory. As is well known to the SD experts, the 
form of the solutions obtained depends crucially 
on the way one chooses to model the QCD run- 
ning coupling at low energies. In this talk we 
summarize results presented in |3] on the impact 
of the analytic invariant charge |3] on CSB. This 
is accomplished by studying the solutions that 
emerge from the standard gap equation governing 
the quark propagator, when the aforementioned 
model for the infrared behavior of the QCD cou- 
pling is adopted. 
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Figure 1. The analytic invariant charge (QJ and 
the perturbative running coupling at the one-loop 
level (solid and dashed curves, respectively). 



We next proceed to the study of the gap equa- 
tion. According to the standard lore, the starting 
point is to express the fully dressed quark propa- 
gator in the following general form [5] 



m 



E(p) = itfA(p 2 ) + B( P 2 ). (2) 



We consider the case without explicit CSB, i.e., 
bare mass mo = 0, and the quark mass is gener- 
ated exclusively through dynamical effects. Then 
CSB takes place when the self-energy S(p) de- 
velops a non-zero value. Traditionally one de- 
fines the quark mass function M(p 2 ) in terms 
of the functions A(p 2 ) and B{p 2 ), as M(p 2 ) = 
B{p 2 ) / A(p 2 ); then CSB occurs when B(p 2 ) ^ 0. 

The quark self-energy, which is represented 
schematically in Figure can be written as 



2(P) = 3^ 



d 4 q 



^ lllS (q)T v (q,p)D^(k), (3) 



where we have used that J2 a ^ a ^ a = V3, 
A a being the Gell-Mann matrices, and k = 
p — q. According to this equation, the self- 
energy is dynamically determined in terms 
of itself, the full gluon propagator, denoted by 
Df v {k) = 5 ab D^(k), and the full quark-gluon 
vertex T u (q,p). Of course, both D°$,(k) and 
T u (q,p) obey their own complicated SD equation, 
a fact which eventually makes unavoidable the use 




Figure 2. The SD equation J2J for the quark 
self-energy. The black blobs represent the fully 
dressed quark and gluon propagators and the 
white one is the proper quark-gluon vertex. 



of simplifications and further modeling of the un- 
known functions involved. 

The first step toward a construction of a man- 
ageable system of equations is to ignore the ghost 
contributions in the quark SD equation. The lat- 
ter are contained in the full gluon propagator 
and the full quark-gluon vertex. More specifi- 
cally, the full quark-gluon vertex should satisfy 
the Slavnov- Taylor identity 

ifcT^fog) [l + b(k 2 )] = 

= [1 - B(k, q)] S-\p) - S-^q) [1 - B(k, q)], (4) 

where b(k 2 ) is the ghost self-energy and B(k,q) 
is the ghost-quark scattering kernel. Without 
ghosts, equation becomes identical to the text- 
book QED Ward identity and the quark-gluon 
vertex now obeys 



ik^T^(p, q ) = S- 1 ip)-S- 1 (q). 



(5) 



This last Ward identity enforces the equality 
Z\ = Z2 between the vertex- and quark-wave- 
function renormalization constants, exactly as 
happens in QED. Evidently, in this treatment we 
are effectively assuming a theory with Abelian- 
like characteristics, where it is hoped that the 
form of the effective gluon propagator employed 
retains some memory of the omitted ghost effects. 

The next step is to employ a "gauge tech- 
nique" inspired Ansatz 7] for the unknown vertex 
T^(p,q). According to this technique, r M (p, q) is 
written in terms of S(p) in such way as to satisfy 
by construction the Ward identity Clearly, 
this procedure introduces the additional ambigu- 
ity of how to fix the transverse part of the vertex, 
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which may lead in higher orders to mishandling 
of the overlapping divergences, but is expected to 
be of little consequence in the infrared. In the fol- 
lowing analysis we will use for r ^ip, q) the simple 
Ansatz proposed in references 6 I 



k 2 



{z[A(p 2 ) - A(q 2 )] ^ + [B(p 2 ) B(q 2 )}}.(6) 



Another important consequence of the forced 
"Abelianization" imposed on the theory estab- 
lishes finally the required link between the gap 
equation and the effective QCD charge. In par- 
ticular, since now Z\ = Z2, one may define a 
renormalization-group-invariant quantity, to be 
denoted by a(k 2 ), which is the exact analogue 
of the QED effective charge, namely 



g 2 D llv {k) = 



k^k v \ 47ra(fc 2 ) 
k 2 J k 2 ' 



(7) 



Choosing the Landau gauge leads to the further 
simplification 



g 2 D^(k)T v (q,p) = g 2 D^(k)A(q 2 ) lv , 



(8) 



since, in that case, the gluon propagator is com- 
pletely transverse. 

Substituting equations and JHJ into the 
quark gap equation J3J), one arrives at the com- 
monly used coupled system for the the quark self- 
energy in terms of A(p 2 ) and B(p 2 ) |5I8| . The last 
necessary step in order to set up our final quark 
SD equation is to employ the usual angular ap- 
proximation, which allows to handle the depen- 
dence on the angle appearing in the arguments of 
the various functions entering into the gap equa- 
tion. Even though the robustness of this approxi- 
mation has been occasionally put into doubt (see, 
for example, |9]), it offers the major advantage of 
reducing the coupled system to one single equa- 
tion, since it automatically forces the relation 
A(p 2 ) = 1. Therefore, we can straightforwardly 
relate B(p 2 ) to the dynamical mass M(p 2 ): 
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(9) 



where M(x) = M(xA 2 )/A, x = p 2 / A 2 , and y = 
q 2 /A 2 . 

Evidently, the invariant charge a(k 2 ) lies in the 
heart of the above equation. Indeed, after the 
succession of simplifications listed above, it is the 
sole ingredient which remains to be modeled. It 
is at this point that the analytic charge ain (A: 2 ) 
of Q will enter into the game, being plugged into 
equation as our model for the QCD effective 
charge a(k 2 ). 

Before solving the resulting integral equation 
numerically, we can infer some qualitative con- 
clusions about the deep infrared behavior of the 
solutions by differentiating both sides twice with 
respect to x (see also references 1618) 1. thus con- 
verting JHJ into a differential equation 



-00,(10) 



where v — lnx. This equation may be solved 
iteratively; the first iteration, corresponding to 
the leading infrared behavior of the dynamical 
mass function M (x) , is obtained by omitting the 
second-order derivative in equation JTDJl. Then, 
the solution is 



Mix) 



In I lnx|, 



0. 



(11) 



The above approximate solution presents a singu- 
lar behavior in the low-energy region, which can 
be interpreted as a hint for confinement (see ad- 
ditional discussion in references |4I8I1()) 1. 

As for the ultraviolet behavior of the solution, 
the conservation of the axial- vector current even- 
tually leads, for sufficient large momenta, to 



Mix) ~ —(In a;) 

x 



A-l 



OO, 



(12) 



where A = 12/ (33 — 2n{) is the anomalous dimen- 
sion of the mass, and D is a constant independent 
of the renormalization point and directly related 
to the quark condensate (see also review jS] and 
references therein). 

The dynamical quark mass function Mix), 
obtained solving numerically the integral equa- 
tion @, is presented in Figure |21 Indeed, one 
can see a soft enhancement of M(x) when x — ► 0, 
as forecasted by Eq. 1 )11)1 . In order to restore a 
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Figure 3. The dimensionless quark dynamical 
mass A4 (x) , obtained from Eq. 10 using the one- 
loop analytic charge JIJ for m = 4 active flavors. 



physical scale in the dimensionless mass function 
A4(x), one has to relate the obtained solution to 
a QCD observable. In particular, this can be ac- 
complished by making use of the method devel- 
oped by Pagels, Stokar JT], and Cornwall [TU|. 
which provides the required link between A4(x) 
and the pion decay constant: 



, 2= 3Al 



2 dy 



yM{y)dy 



[y + M*(y)Y 



.(13) 



Thus, the scale parameter A (which is the only ad- 
justable parameter within the approach at hand) 
can be evaluated by requiring the right-hand side 
of Eq. (|13(l to acquire the experimental value of 
the pion decay constant fl = (93MeV) 2 . Ul- 
timately, for the case of m = 4 active flavors, 
this results in estimate A = 880 MeV. The higher 
loop corrections to the analytic running coupling 
Q an (/: 2 ) [Hj do not alter qualitatively the picture 
obtained. Specifically, the confining behavior of 
the dynamical mass function A4(x) persists, but 
the infrared singularity displayed is weaker than 
that of Eq. . In turn, this further elevates the 
estimate of the scale parameter A. 

The conclusion that we draw from the above 
analysis is that the analytic invariant charge de- 
veloped in Ref. [3] and the SD equations may co- 
exist in a complementary and qualitatively con- 



sistent picture, and that the obtained value for 
the only free parameter turns out to be in the 
right ballpark. 
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